Reentrant Kondo effect in Landau quantized graphene 
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We have studied the interplay of an Anderson impurity in Landau quantized graphene, with special 
■ emphasis on the influence of the chemical potential. Within the slave-boson mean-field theory, we 

' found reentrant Kondo behaviour by varying the chemical potential or gate voltage. Between Landau 

levels, the density of states is suppressed, and by changing the graphene's Fermi energy, we cross from 
£SJ i metallic to semiconducting regions. Hence, the corresponding Kondo behaviour is also influenced. 

The f-level spectral function reveals both the presence of Landau levels in the conduction band and 
the Kondo resonance. 
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I. INTRODUCTION 

qj • The existence of Dirac fermionic excitations in single layer graphene has attracted great interes!^^^ Theoret- 
ically the linear dispersion and chiral nature of massless quasiparticles has many unexpected physical consequences, 
especially for magnetotransport properties'^ ' 10 ' 11 . Experimentally graphene may be considered as a zero gap semi- 
conductor with charge carriers of very high mobility whose density may easily be controlled by applying a gate voltage 
"j^ . which gives hope for interesting applications 6 -. 

The scattering of graphene quasiparticles from normal impurities determines electronic and thermal transport and 
I , is well studied ' 10 ' 11 ! 12 ' 13 . The linear dispersion of quasiparticles and the associated pseudogap plays an essential 
role in the magnetotransportii because it leads to unconventional Landau level quantization as compared to the case 
of parabolic bands. One essential difference is that a zeroth order Landau level is pinned at zero energy for any field 
strength. With increasing impurity scattering the pseudogap is therefore gradually filled and the density of states 
exhibits oscillations as function of the chemical potential or carrier number^ii. 

Because the graphene carrier number and density of states can easily be manipulated with gate voltage and magnetic 
field it is attractive to consider the effect of magnetic impurities in a graphene host. It is well known that magnetic 
moments in the Fermi sea show the competition of on-site Kondo singlet formation and intersite RKKY interactions 
^ji^ ' which is highly sensitive to the density of states around the Fermi level^. While experiments in graphene are still 
, missing we think that chemisorption or adsorption of transition metal ions on graphene sheets is the most likely 
—i ■ realization of Kondo effect and RKKY interaction with Dirac electrons. Possible ways for creation of local moments 

, ' in carbon based materials have been discussed in Refs. mmn 

■ Theoretically the Kondo and Anderson- type models for local moments in graphene have been studied without 
' magnetic fiel d 16 ' 18 . In the present work we investigate the effect of Landau level formation on the screening of 

f^) , magnetic impurities described by the Anderson model. Due to the strong DOS variation with chemical potential 

^ ' and magnetic field it should be possible to drive the magnetic impurity in and out of the Kondo regime which is 

• i-H , characterized by the formation of an f-level resonance close to the chemical potential. This effect is studied within 

■ mean field slave boson theory and assuming a constant broadening of Landau levels by normal impurities. In Sect. II 
J-j | we define the Anderson impurity model coupled to a Dirac Fermi sea in a magnetic field. The solution of the model 

. - - i within the saddle point approximation is derived in Sect. Ill and the numerical results are discussed in Sect. IV. 
Finally Sect. V gives the conclusion. 
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II. ANDERSON IMPURITY IN GRAPHENE 

The Hamiltonian of Dirac quasiparticles living on a single graphene sheet, interacting with an infinite-?/ Anderson 
impurity at the origin is given b y 7 ' 19 i 20 : 

a 

(i) 

where Cj's are the Pauli matrices, and stand for Bloch states residing on the two different sublattices of the bipartite 
hexagonal lattice of graphenoli^. The quasiparticle spectrum vanishes at six point in the Brillouin zone. Out of these 
six, only two are non-equivalent points, and arc referred to as K and K' points in the Brillouin zone, denoted by the s = 
1 and -1 index, respectively, <5 is the Kronecker delta. The Kronecker delta function accounts for the non-equivalence 
of the two Dirac points^ at K and K' . The finite chemical potential /i accounts for particle- hole symmetry breaking. 
The vector potential for a constant, arbitrarily oriented magnetic field reads as A(r) = (—By cos 9, 0, B(y sin 9 cos </> — 
x sin 9 sin </>) ) , where 9 is the angle the magnetic field makes from the z axis, and </> is the in-plane polar-angle measured 
from the x-axis. The Zeeman term is assumed to couple to both the impurity and Dirac electrons by the same g 
factor, h = gpiBB, vp w 10 6 m/s, is characteristic to graphene. /+ and f a creates and annihilates an electron on 
the localized E level, b + and b are the slave boson operators, responsible for the hole state o 21 i 22 ' 23 . These take the 
infinite- U term into account. 

The Hamiltonian should be restricted to the subspace 

E rtf«+ b+b = L ( 2 ) 

a 

Within the mean-field approximation, the slave-boson operators are replaced by their expectation value, bo = (b), 
and the constraint is satisfied by introducing a Lagrange multiplier A: 

# CT)S (r) + E(£ + A-/ l a)/+/ CT + 
+Vb Q £ (/+*„,.((>) + *+ s (0)/ CT ) + A(6 2 - 1). (3) 

£T,S 

In the particle-hole symmetric case (fi = 0), in the absence of magnetic field, the energy spectrum of the system of 
Dirac electrons is given by 

E(k)=±v F \k\. (4) 

This describes masslcss rclativistic fcrmions with spectrum consisting of two cones, touching each other at the end- 
points. From this, the density of states per spin follows as 

pH = -^%-B(k)) = ^ [kdkS(u>±vpk) = ^-, (5) 

7T k 7T IT J 

where k c is the cutoff, D — vpk c is the bandwidth, and A c — 47r/fc 2 is the area of the hexagonal unit cell. We mention 
in passing, that an applied gate voltage directly controls the number of extra charge carriers in the system, which is 
given by eV ~ n = p(uj)dtjj = fi 2 / D 2 . Hence, chemical potential is proportional to the square root of the applied 
gate voltage even in the ideal case, without any scatterers and magnetic field. Such a relation can hardly be calculated 
for the realistic case. Nevertheless, the chemical potential always varies monotonically with the gate voltage due to 
the positiveness of the density of states. 

Magnetic impurities in gapless fermi systems have thoroughly been studied starting with the pioneering work of 
Withoff and Fradkin^i, and the focus was on the influence of gapless excitations on the Kondo phenomenon (for 
a review, see Ref. [25h . Recently, Kondo effect in graphene without magnetic field has been studied within this 
frameworki£. Here we allow for Landau quantization of the quasiparticle spectrum, and study the orbital effect of 
magnetic field on the various Kondo phases. 



H = j dr£*+,(r) 



— vp(l — 28jy8 s _i)<Tj {—idj + eAj(r)) — /.t — ha 
(/+* CTlS (0)6 + &+*+ s (0)/ CT ) , 



H = 



(r) 



—Vf(1 - 2dj, y 5 s ,-i)<7j (-idj + eAj(r)) - fi - ha 
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In the presence of magnetic field, the eigenvalue problem of our Hamiltonian without the localized level can readily 
be solved^. From now on, we concentrate on the K point, the eigenfunctions of the K' point can be obtained by 
exchanging the two components of the spinor. Momentarily, we also neglect the spin and the Zeeman term, and 
concentrate on the effect of Landau quantization. They will be reinserted when necessary. For the zero energy mode 
(E=0), the eigenfunction is obtained as 

e ikx / q 



and the two components of the spinor describe the two bands. The energy of the other modes reads as 



E(n, a) = auj c \/ n + 1 (7) 



with a = ±1, n = 0, 1, 2,.. . , u> c = vpy / 2e\B cos(#)| is the Landau scale or energy, but is different from the cyclotron 
frequency 2 ^. Only the perpendicular component of the field enters into these expressions, and by tilting the field away 
from the perpendicular direction corresponds to a smaller effective field. The sum over integer n's is cut off at N 
given by N + 1 = (D/uj c ) 2 , which means that we consider 2N + 3 Landau levels altogether. 
The corresponding wave function is 



with cyclotron length It, = 1/V eB. Here 4> n {x) is the nth eigenfunction of the usual one-dimensional harmonic 
oscillator. The electron- field operator at the K point can be built up from these functions as 



,a.C-k,n,a 



(9) 



The Green's functions of these new operators do not depend on fc, and read as 



G (lUJ n ,k) = — , (10) 
lU) n 

G (iw n , k, n, a) = l — (11) 

iuj n — tL[n^ a) 

for Cfc and Ck, n ,a, respectively, and oj n is the fermionic Matsubara frequency. As seen from above, the density of 
states in the presence of quantizing magnetic field contains Dirac-delta peaks located at the Landau level energies, 
and quasiparticle excitations have infinite lifetime. To describe a more realistic situation, scattering from disorder 
needs to be considered in the presence of the magnetic field. Usually, the resulting self energy of the Dirac fermions, 
determined in a self-consistent manner, depends on the frequency and field strength^ii. However, good agreement 
can be reached by assuming a constant, phcnomcnological scattering rate, denoted by T, for small and moderate 
fields, as can be learned from similar analyzi o 9 ' 27 . To simplify calculations, we have chosen to mimic disorder by a 
constant scattering rate. 



III. SADDLE-POINT EQUATIONS 



The free energy of the system can be found from Eq. using standard technique^. The value of A and bo is 
determined self-consistently by minimizing the free energy of the system with respect to the m 28 i 29 . As a result, by 
restoring the finite chemical potential, spin and the Zeeman term, the saddle point equations at T = are given by: 

b o= Si{E + \-ah^ + ah), (12) 

(T = ±l 

A= S 2 (E + \-ah,fj, + ah), (13) 

cr=±l 



where the auxiliary functions arc defined as 

oo 

dx 



Si(y,z) 



y + &§2V 2 ImE(a; + T 



S 2 (y,z) = 2V 2 



E(*) = 



o 

f,2oT/2 



it (x + bl2V 2 ReT,(x + T- iz)) 2 + (y + & 2 2^ 2 Im£(.T + T — iz)) 2 ' 



(14) 



dx (x + b^2V 2 ReT,(x + T- iz))ReE(x + T - iz) + ImE(x + T — iz)(y + b 2 2V 2 lmT,(x + T- iz)) 



(x + bl2V 2 ReY,(x + T — iz)) 2 + (y + 6§21/ 2 ImS(x + T — iz)) 2 



1 



7V + 1 



N 

E 

fc=0 



2: 



Z 2 +UJ 2 {k+l) 



1 



7V + 1 



iV + 2 



(15) 
(16) 



The 1/(N + 1) oc B prefactor denotes the Landau level degeneracy. These reduce to the standard saddle-point 
equations for gapless phase o 28 ' 30 for zero field and T = 0. The extra factor of 2 in front of V 2 stems from the two 
non-equivalent Dirac cones at the K and K' points, since by using 



dk<f> n (y - kl 2 B )(j> m (y - kl 2 B ) 



I 2 

' B 



(17) 



each Landau level at each valley contributes to the hybridization energy by V 2 . Eq. (j 12[) accounts for the constraint 
of having at most one f-electron at the impurity site (Eq. ([2])), Eq. (fT3|) stems from the equation of motion of the 
slave boson field b: since it is constant in the mean-field approach, the terms determining its dynamics should add 
up to zero^ 2 -. Here S(z), which is related to the f-level self energy, contains all the information about the conduction 
electron bath, where the magnetic impurity is embedded. When the field strength goes to zero, these equations reduce 
to those found in gapless phases 2 ^ as 



= 2 ^>(l + ^ r 



(18) 



In this case, for fi = V = 0, the critical f-level energy is found to be E c = —8V 2 /D. For E < E C1 only the trivial 
solution occurs (bo = 0), hence charge fluctuations are completely suppressed. The solution of Eqs. (fT2|) -([T3 ll can be 
carried out by realizing, that Eq. (|12[) depends only on the renormalized f-level energy E + A, and not separately on 
the two variables E and A. Then, by fixing the value of E + A, we can determine the corresponding bo by iteration, 
for example. By inserting the values of the renormalized f-level energy and the slave-boson expectation value to Eq. 
(|13[) . we can directly read off A, and determine E by subtracting it from the renormalized f- level energy. As in other 
similar problems, this method predict a quantum phase transition at T = 0, characteristic to gapless Kondo phases. 
However, in our case, the order of the transition can change from second to first. Such a transition is absent for 
magnetic impurities embedded to normal metals 2 ^. 



IV. SPECTRAL FUNCTION, DISCUSSION 



The full f-clectron spectral function per spin along the real frequency axis reads as 



Pf,a{u) = Im — 

7T UJ 



It 



E -X + ah- 2V 2 b 2 Y, f {uj + fx + ah + iT) ' 



(19) 



where S/(a;) = — ix). It shows the Landau level oscillations, and in addition, a large Kondo peak whenever 
non-trivial solution of Eqs. (|12p - (|13p exists. Without magnetic field and normal impurities, it shows marginal Fermi 
liquid behaviour—, as can be observed from the analytically continued f-level self-energy: 



4^ 2 6 2 
D 2 



( — 



(20) 



which, in the presence of normal impurities, turns into the usual fermionic self energy at low energies, with u> 
being replaced by T, the scattering rate. In the case of quantizing magnetic field, without normal impurities, its 
imaginary part consists of Dirac-delta functions at the Landau level energies. Non-magnetic impurities smear these 
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FIG. 1: (Color online) The order parameter (left panel) and the contour plot (for _E/D=-0.036, right panel) of the f-level 
density of states are shown for N = 1000, 2(V/D) 2 = 0.01, T — 0.05o; c , h — 0.1uj c . For negative chemical potentials, the order 
of the transition changes from second to first, as is denoted by the thick black line. The parallel structures in the contour plot 
denote the Kondo peaks, separated by twice the Zeeman energy. 

singularities, by transforming them into Lorentz functions. Hence, the self energy is that of a well-formed Fermi liquid 
for realistic situations. We mention the possibility of detecting non-Fermi liquid phase a 15 i 16 in graphene due to the 
valley degeneracy, which might originate from multichannel Kondo physics. This can be signaled by the finite critical 
value of E even at T — as is seen in Fig. [T] However, such a situation is unlikely to be observed by our simple 
mean-field analysis. 

The numerical solution of the saddle-point equations have been carried out, and the result for the order parameter 
(&o) is shown in Fig. Q] It is directly related to the f-level occupation through Eq. ((2|) as nt = 1 — 6g. In general, 
by increasing a Fermi surface develops, and the critical value of E penetrates deeply into the E < regional. 
However, there is a crucial difference between positive and negative chemical potentials: for the former, the Kondo 
energy scale becomes extremely small, as can be seen in Fig. [T]in the f-level density of states, since the f-level energy 
is well below the Fermi energy of conduction electrons. In the latter case, the energy level of f-electrons lies above the 
Fermi energy, hence the Kondo scale enhances. 

When the chemical potential is close to a Landau level energy E(n, a), the non-trivial solution (nonzero bo) extends 
further in the E < region, because the density of states is enhanced around Landau level energies. Between 
Landau levels, the conduction electron density of states is suppressed, and the critical f-level energy moves closer 
to zero, leading to oscillatory behaviour in the phase boundary and diagram, as can be seen in Fig. [TJ Here we 
assume N = 1000, which corresponds to weak or moderate fields^, depending on the explicit value of the cutoff 
D. The Zeeman term is chosen to be much smaller than u) c , is follows from actual numbers in graphene^. As a 
result, by changing the chemical potential, we can move between Landau levels, and we can enter into and leave the 
Kondo regime. Hence reentrant behaviour is found. By decreasing the scattering rate T from normal impurities, the 
oscillations along the phase boundary become more pronounced. 

The reentrant behaviour can more directly be checked in the f-level density of states in Fig. [2l where large Kondo 
peaks are observable (due to Zeeman splitting) in addition to small oscillation stemming from Landau levels in the 
conduction band, when non-trivial solution exists. The distance between the two parallel ridges in the contour plot 
in Fig. [T]is two times the Zeeman energy, as it should be. When only the trivial solution exists (bo = 0), Pf(u>) is 
completely suppressed, since the maximally allowed one particle always occupies the f-level. 

In the Kondo regime, as one varies the chemical potential, the Kondo temperature docs not change monotonically. 
It remains mainly pinned to the closest Landau level, and then suddenly jumps to the neighbouring one, as is seen in 

Fig. m 

These features in the density of states can probably be detected by conductance measurements, which measures 
directly the inverse of the f-level density of states, in addition to normal impurities — . When pf(0) is finite, its 
contribution is thought to overwhelm that of normal impurities^. Then, both the Landau level oscillation stemming 
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FIG. 2: (Color online) In the top left panel, the spin averaged f-level density of states is shown as a function of the chemical 
potential for TV = 1000, 2(V/D) 2 = 0.01, T = 0.05o; c , h = 0.1w c for £/£>=-0.036 (blue dashed), -0.0387 (red solid) and -0.04 
(black dashed dotted). Note the presence of states for the middle value close to zero chemical potential. Reentrant behaviour 
is also observable close to the first Landau level for _E/D=-0.04, fj, ~ ui c - The three 3 dimensional plots show the evolution of 
the spectral density as a function of chemical potential and frequency for the previous three E values from top to bottom, left 
to right. Note the presence of small island of states for E/D = —0.0387, responsible for the Kondo effect at fj, « 0. 



from the orbital quantization of conduction electrons and the Kondo behaviour could be seen. The change in the 
Kondo temperature (the peak position on Fig. [T]and[5J) by chemical potential or gate voltage can also in principle be 
detected from the specific heat. This is expected to exhibit a double peak structure around Tr due to the Zeeman 
term. However, such measurement on thin graphenc films are extremely demanding. The change of the spectral 
function would reveal itself directly in photoemission spectroscopy^, which, however, in a finite magnetic field, does 
not constitute a standard experiment. 

The presence of Kondo resonance makes itself felt in magnetic responses, which are expected to be different from the 
usual Kondo behaviour due to the presence of orbital (or diamagnetic) terms. This already influences the magnetic 
susceptibility of pure graphene, and probably alters the magnetic behaviour of the impurity as well. 
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V. CONCLUSIONS 



In summary, we have studied the infinite-U Anderson impurity embedded to a host of two-dimensional Dirac 
fermions within the self-consistent slave boson mean field theory. The host material corresponds to graphene, where 
the elementary excitations on a honeycomb lattice are Dirac fermions. Such a system can most probably be realized by 
the chemisorption of transition metal ions on graphene sheets. We allow for the Landau quantization of the conduction 
electron spectrum, which turns out to be unusual^^ (E oc y/n) in comparison with normal metals (£ a n). Then, 
we study the effect of orbitally quantizing magnetic field on the Kondo phenomenon, in addition to the Zeeman term, 
albeit the latter is thought to be negligible in the presence of Landau levels^. 

When the chemical potential lies close to a Landau level energy, the conduction electron density of states is enhanced, 
and the mixed-valence regime extends further in the E < regional. Between Landau levels, the density of states 
resembles to that in an insulator, and the local moment regime gains ground. Hence by varying the chemical potential 
between Landau levels, reentrant Kondo behaviour is found. This manifests itself in the f-electron spectral function, 
which accommodates small islands of states, corresponding to the mixed-valence case, separated by the deserts of 
local moment region, as \i changes. The reentrant behaviour manifest itself strongly around the zeroth order Landau 
level, but for clean samples, it should be observable at higher levels as well. 

The experimental detection of this phenomenon can be done by conductance measurements in magnetic field at 
low temperatures, which can show the presence or absence of Kondo peaks in the f-electron density of states. 
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